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INFLUENCE  OF  COMBUSTION  PROCESS  ON  STABILITY* 
by  Richard  J.  Priem 
Lewis  Research  Center 

SUMMARY 

I  The  relative  importance  of  various  steps  in  the  combustion  process  on  combustion 
instability  was  compared  theoretically  by  assuming  that  individual  steps  controlled  the 
growth  of  a  disturbance.  The  combustion  steps  involved  the  time  dependence  of  (1)  pro¬ 
pellant  injection,  (2)  accumulation  of  unburned  propellant,  (3)  propellant  atomization, 

(4)  propellant  vaporization,  and  (5)  chemical  reaction.  The  calculation  gave  stability 
limits  (minimum  disturbance  that  will  grow  to  instability)  with  various  steps  considered 
to  control  the  combustion  rate. 

The  calculations  show  that  chemical  reaction  is  the  most  sensitive  mechanism,  but 
if  accumulation  of  unburned  propellants  is  included,  the  chemical  reaction  step  is  likely 
to  be  important  only  for  low  conversion  rates  that  do  not  occur  in  most  rocket  combus¬ 
tors.  With  realistic  combustion  conditions  the  physical  processes  of  vaporization  and 
atomization  appear  to  be  the  most  important  in  determining  instability  limits,  the  limits 
for  each  step  depending  on  the  accumulation  of  unburned  propellant  and  vortex  velocity. 
Experimental  instability  studies  with  an  8-foot  torus  combustor,  which  is  similar  in 
geometry  to  that  used  for  the  theoretical  calculations,  confirmed  the  calculated  results. 
Experiments  also  confirmed  qualitatively  the  fluctuations  that  occur  during  instability  J 

INTRODUCTION 

Combustion  of  propellants  in  a  liquid  rocket  combustion  chamber  involves  the  re¬ 
action  of  a  reducing  agent  (fuel)  and  an  oxidizing  agent  to  produce  reaction  products  with 
the  liberation  of  heat.  This  seemingly  simple  statement  of  a  reaction  between  the  reduc¬ 
ing  and  oxidizing  agents  can  become  very  complicated  when  the  mechanism  or  steps  are 
considered  that  occur  between  the  injection  of  liquid  propellants  and  the  time  or  position 
that  the  final  products  are  formed.  Occurring  simultaneously  with  the  combustion  are 
fluctuations  in  the  chamber  pressure,  which,  under  certain  conditions,  are  seen  to 
amplify  into  periodic  oscillations  driven  by  the  combustion  process.  These  periodic 
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Figure  1,  -  Combustion  section  considered  for 
various  instability  models. 


pressure  oscillations  result  in  engine  failure  due  to  ex¬ 
cessive  heat  transfer  to  the  walls,  which  eventually 
leads  to  thermal  failure,  or  vibrational  interference  with 
controls  and  other  accessories  that  leads  to  mechanical 
failure.  Thus,  it  is  important  to  understand  the  role  of 
the  combustion  process  in  the  growth  of  these  oscilla¬ 
tions. 

Previous  studies  have  shown  how  various  parts  of 
the  complicated  combustion  process  could  support  the 
growth  of  a  pressure  disturbance.  (See  bibiography  in  ref.  1. )  These  studies  have 
usually  resulted  in  qualitative  descriptions  of  the  combustion  process  with  pressure 
oscillations  or  a  theory  for  a  particular  step  that  could  sustain  an  oscillation.  In  each 
study  the  assumptions,  combustor  design,  model,  etc. ,  are  different;  thus,  it  is  im¬ 
possible  to  evaluate  the  importance  of  a  particular  step  in  the  combustion  process  in  sup¬ 
porting  the  growth  of  a  pressure  disturbance.  Therefore,  a  theoretical  study  was  con¬ 
ducted,  similar  to  that  reported  in  reference  2  which  used  a  toroidal  ring  to  represent 
the  combustor,  in  which  the  various  combustion  steps  were  individually  assumed  to  con¬ 
trol  the  growth  of  a  disturbance.  It  was  determined  which  step  could  amplify  the  small¬ 
est  disturbance,  so  that  this  process  could  be  further  studied  under  more  realistic  and 
complicated  conditions.  To  determine  the  influence  of  a  vortex  in  the  combustor,  which 
may  be  produced  by  viscous  forces  with  a  transverse  oscillation  (ref.  3),  calculations 
were  also  made  with  a  constant  vortex  velocity  superimposed  on  the  oscillations.  As  a 
check  on  the  theory,  the  calculated  stability  limits  were  compared  to  the  experimental 
limits  obtained  with  several  injectors  in  a  torus  combustor. 


THEORY 

General  Equations 

The  combustor  geometry  and  transport  equations  used  in  the  analytical  portion  of 
this  investigation  are  identical  to  those  in  reference  2.  The  combustor  was  an  annular 
section  with  a  very  small  thickness  AR  and  length  Az,  as  shown  in  figure  1.  This 
geometry  limits  the  oscillation  to  a  transverse  wave.  The  propellants  are  assumed  to 
be  uniformly  injected  and  to  burn  at  some  rate  which  depends  on  their  position  within  the 
combustor.  The  burning  is  accompanied  by  gas  and  liquid  flow  in  the  axial  direction. 

As  the  propellants  burn,  random  disturbances  occur  within  the  combustor.  If  the  ampli¬ 
tude  of  a  disturbance  is  large  enough,  it  may  develop  into  a  wave  in  the  combustor .  The 
object  of  this  study  was  to  determine  analytically,  for  various  combustion  steps,  the 
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minimum  amplitude  of  a  pressure  disturbance  in  the  annulus  that  is  needed  to  develop 
a  wave.  This  calculation  was  accomplished  by  solving  the  transport  equations  that 
follow,  which  are  derived  in  reference  2.  Cylindrical  coordinates  are  used  in  which  the 
z -direction  corresponds  to  the  axial  direction  in  the  combustor. 
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(Symbols  are  defined  in  appendix  A. )  The  derivatives  in  the  axial  (z)  direction  were 
determined  with  the  assumption  that  the  total  mass,  momentum,  and  energy  in  the 
annulus  were  constant.  Furthermore,  it  was  assumed  that  these  derivatives  were  in¬ 
dependent  of  0.  These  assumptions  result  in  the  following  equations  (derived  in  ref.  2), 
which  permit  the  evaluation  of  the  derivatives  with  respect  to  z: 
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Combustion  Models  for  Burning  Rate 

In  addition  to  the  transport  equations,  an  expression  for  the  instantaneous  local 
burning  rate  is  required.  A  schematic  model  for  the  entire  combustion  process  is 
illustrated  in  figure  2.  Propellants  are  assumed  to  enter  the  chamber  at  a  varying  rate 
co^  dependent  on  the  injector  pressure  drop  -  P^,  which  assumes  no  inertial  or 
capacitance  factors  in  the  feed  system.  (The  derivations  of  expressions  for  various 
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Figure  2.  -  Schematic  model  of  combustion  process. 


rates  and  masses  are  given  in  appendix  B. )  Liquid  is  atomized  at  a  rate  dependent 

a 

on  th^  amount  of  unatomized  liquid  M  .  Similarly,  the  propellants  vaporize  at  a  rate 

a 

0)^  depending  on  the  amount  atomized  but  not  vaporized  and  then  react  chemically 
at  a  rate  depending  on  the  amount  of  each  propellant  vaporized  but  unreacted  M^. 

In  figure  2  the  arrows  represent  propellants  changing  from  one  state  to  another,  and  the 
crosshatched  areas  indicate  propellants  in  a  given  state.  The  steps  between  propellant 
vaporization  and  reaction  must  include  mixing,  but  since  there  are  no  known  expressions 
to  determine  the  local  instantaneous  mixing  rate,  it  was  not  included  in  the  model. 

To  determine  the  relative  importance  of  the  various  combustion  steps  in  driving  a 
pressure  oscillation,  several  models  were  used.  Each  model  assumed  that  one  of  the 
mechanisms  controlled  the  local  instantaneous  combustion  rate  and  that  the  other  mecha¬ 
nisms  exerted  negligible  effect  on  the  overall  rate;  subsequent  mechanisms  were  con¬ 
sidered  very  fast  so  that  they  would  follow  at  the  same  rate.  Preparatory  mechanisms 
were  considered  fast  enough  to  follow  the  injection  rate  and  yet  be  insensitive  to  local 
fluctuations.  With  each  model  an  adiabatic  pressure  disturbance  was  applied  to  the  sys¬ 
tem,  and  its  growth,  or  decay,  was  observed.  This  disturbance  was  given  by 


=  1  +  sin  0 


(y-l)/y 


T'  =  1  +-P  sin  6’ 

V  2  y 


(9a) 


(9b) 


5 


/  p  \ 

p»  =  h  +  ^  sin  e’j 


(9c) 


Chemical  Reaction  Model 

For  this  model  it  was  assumed  that  only  the  chemical  reaction  rate  varied  during  in 
stability.  The  burning  rate  expression  was  then 
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Vaporization  Model 

For  this  model  it  was  assumed  that  the  vaporization  rates  for  oxidizer  and  fuel  were 
equal  and  were  the  slowest  compared  to  those  of  the  other  rates.  It  was  also  assumed 
that  the  atomization  and  chemical  reaction  mechanisms  were  very  rapid  so  that  they 
followed  the  injection  and  vaporization  rates,  respectively,  or 


and 
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The  burning  rate  expression  then  becomes 
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Atomization  Model 

It  was  assumed  that  the  atoraisalion  rates  for  oxidizer  and  fuel  were  equal  and  were 
the  slowest  for  this  model  and  that  the  vaporization  and  chemical  reaction  followed  them. 
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The  burning  rate  was  then  obtained  from 
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Burning  rate  parameter,  2? 

Figure  3.  -  Instability  boundaries  for  chemical 
reaction  model  with  finite  concentrations. 


THEORETICAL  RESULTS 

The  results  of  the  numerical  calculations  were  pressure,  particle  velocity,  density, 
temperature,  burning  rate,  and  concentration  of  mass  at  each  of  20  positions  in  the 
annulus  as  functions  of  time.  Typical  plots  of  pressure  and  velocity  can  be  found  in 
reference  2.  From  these  types  of  curves  it  was  possible  to  determine  the  amplitude  of 
the  minimum  disturbance  that  would  increase  with  time.  This  value  of  the  minimum  dis¬ 
turbance  was  used  to  determine  the  boxmdaries  between  stable  and  unstable  combustion. 
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All  calculations  presented  herein  were  for  the  following  conditions  of  v^,  AP', 
Av',  and  7  in  equations  (1)  to  (17); 

Viscous-dissipation  parameter,^ . 

Nondimensional  axial  gas  velocity,  v’ . 

La 

Nondimensional  steady-state  pressure  drop  across  injector,  AP' 

Nondimensional  axial  velocity  difference,  Av' . 

Specific-heat  ratio,  y . 

The  influence  of  these  parameters  on  the  calculated  results  can  be  found  in  reference  2. 


0. 

.  0.05 

.  .  0.2 
.  0.01 
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Chemical  Reaction  Model 

stability  limits  for  the  chemical  reaction- rate  model  are  shown  in  figure  3  for  a 
limited  supply.  With  an  unlimited  supply  all  disturbances  that  were  tested  excited  in¬ 
stability  under  all  operating  conditions  values)  because  of  the  extreme  sensitivity  of 
the  exponential  function  of  the  chemical  reaction  equation  to  a  disturbance.  However, 
when  a  realistic  value  for  the  concentration  of  unburned  propellant  is  used  (0. 1  or  0.  01 
mass  fraction),  the  stability  limits  occur  at  very  low  burning  rate  parameters  =  0. 1 
or  less).  As  discussed  in  reference  2,  these  low  burning  rate  parameters  occur  only 
for  specialized  small-scale  research  combustors. 


Vaporization  Model 

stability  boundaries  for  the  vaporization  model  are  shown  in  figure  4.  The  boundary 
for  infinite  supply  (M,,  ^/co.  ^t  =  a)  is  the  same  as  that  shown  in  reference  2  for  a 
value  of  Av'  of  0.  01.  As  the  quantity  of  unvaporized  material  decreases,  the  stability 
limits  are  increased.  (A  larger  amplitude  of  disturbance  Pp  is  required  to  excite  in¬ 
stability.  )  Since  M  /w.  t  is  equivalent  to  the  ratio  of  the  time  to  burn  an  average 
particle  of  propellant  to  the  wave  time,  it  is  seen  that,  when  the  burning  time  is  much 
less  than  the  wave  time,  for  example,  o^*^i  oV  ~  extremely  large  distur  ¬ 

bance  is  required  to  excite  the  instability. 

The  explanation  for  the  increased  stability  with  lower  values  of  ^t^  can 

be  obtained  from  figure  5,  which  shows  the  time  histories  for  the  local  burning  rate,  the 
unvaporized  mass,  and  the  pressure  at  two  values  of  o^^i  oV'  large  values 

of  o^^i  o^w  unvaporized  mass  remains  constant  during  the  disturbance,  while 
for  low  values  of  M  _t  the  unvaporized  mass  decreases  as  the  pressure  and 

V,  O  o  w 
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Vaporization  mass 


burning  rate  increase.  Since  the  burning  rate  as  given  by  equation  (12)  is  proportional 
to  the  square  root  of  the  unvaporized  mass,  the  burning  rate  is  reduced  by  the 
unvaporized- mass  term  during  the  first  1/4  cycle  of  the  disturbance.  This  reduction 
results  in  less  energy  being  added  to  the  wave;  hence,  the  disturbance  decays. 


Atomization  Model 

stability  boundaries  for  the  atomization  model  are  shown  in  figure  6.  Stability 

lines  are  again  shown  for  various  values  of  M  /o).  t  .  Vv'hen  the  stability  boundaries 

a,  o  1 ,  o  w 

for  the  vaporization  model  (fig.  4)  are  compared  with  those  for  the  atomization  model 
(fig.  6),  it  is  seen  that  in  some  regions  the  vaporization  process  is  more  sensitive  to  a 
disturbance  and  that  in  other  regions  the  atomization  process  has  the  greater  sensitivity. 
Therefore,  both  processes  must  be  considered  important  in  controlling  combustion  in¬ 
stability. 

Time  histories  of  unatomized  liquid,  burning  rate,  and  pressure  are  shown  for  the 
atomization  model  in  figure  7.  The  histories  are  similar  to  those  in  figure  6  for  the 
vaporization  model;  that  is,  for  large  values  of  M  system  is  unstable, 
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(b)  Vaporization  mass  accumulation  factor,  1. 

Figure  5.  -  Time  history  of  disturbance  for  vaporization  model.  Burning-rate 
parameter,  1;  nondimensional  peak-to-peak  pressure  disturbance,  0.01. 

and  the  unatoraized  liquid  remains  constant.  For  small  values  of  o^^i  o^w 
system  is  stable,  and  the  unatomized  liquid  decreases  as  the  burning  rate  and  pressure 
are  increased.  This  decrease  in  unatomized  liquid  reduces  the  burning  rate,  which  re¬ 
sults  in  insufficient  energy  to  drive  the  disturbance. 


Vortex  Flow  Model 

Recent  studies  (ref.  3)  have  shown  that  with  viscous  forces  a  transverse  pressure 
oscillation  can  generate  a  vortex  within  the  combustor.  Therefore  to  determine  the 
significance  of  such  a  flow  pattern,  additional  calculations  were  performed  to  determine 
the  influence  of  a  vortex  on  stability.  These  calculations  were  performed  by  using  the 
vaporization  model  with  a  constant  vortex  velocity  v^^  added  to  the  calculated  values  of 
gas  velocity  in  the  tangential  direction,  v^.  Results  of  these  calculations  are  illus¬ 
trated  in  figure  8.  The  stability  limits  are  dramatically  changed  by  a  small  vortex  vel¬ 
ocity.  For  a  given  vortex  velocity  the  system  has  the  same  stability  at  low  burning  rate 
parameters  that  it  had  for  the  vaporization  model.  At  some  critical  burning  rate  param- 


11 


Figure  6.  -  Instability  boundaries  for  atomization  model. 

eters,  dependent  on  the  vortex  velocity,  the  vortex  model  shows  that  almost  any  distur- 
bance  will  grow.  Larger  critical  burning  rate  parameters  are  required  for  smaller  vor¬ 
tex  velocities.  The  effectiveness  of  baffles  can  be  partially  explained  by  the  vortex 
model  results.  The  vortex  velocity  can  be  eliminated  by  using  baffles,  and  the  stability 
of  the  vaporization  model  boundaries  can  thus  be  increased. 

The  vortex  model  can  also  explain  the  results  reported  for  the  investigation  of  refer¬ 
ence  4,  in  which  a  tangential  gas  flow  was  used  to  excite  instability.  The  quantity  of 
tangential  gas  flow  required  to  induce  instability  was  dependent  on  the  operating  condi¬ 
tions  of  the  combustor. 


EXPERIMENTAL  RESULTS 

To  test  the  validity  of  the  theoretical  results,  an  experimental  program  was  also 
conducted.  Since  the  calculations  were  only  performed  for  an  annular  segment  of  a  com¬ 
bustor,  an  annular  or  torus  combustor  was  used  in  the  experimental  program  (fig.  9, 
p.  15).  The  torus  had  an  8-foot  diameter  and  10  injectors  with  nozzles  equally  spaced 
around  the  chamber.  Four  different  types  of  injectors  were  used  to  change  the  operating 
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Figure  7.  -  Time  history  for  atomization  model.  Burning-rate  parameter,  1;  non- 
dimensional  peak-to-peak  pressure  disturbance,  0.02. 


conditions.  The  injectors  (shown  in  fig.  9)  were  (1)  an  unlike-triplet  injector  with  two 
hydrogen  streams  impinging  on  one  liquid-oxygen  jet,  (2)  an  impinging-jet  injector  with 
two  liquid-oxygen  jets  impinging  on  each  other  and  four  parallel  jets  of  hydrogen,  (3)  a 
coarse  parallel-jet  injector  with  one  jet  of  liquid  oxygen  and  two  jets  of  hydrogen,  (4)  a 
fine  parallel-jet  injector  with  four  jets  of  liquid  oxygen  and  four  of  hydrogen.  Two-inch- 
diameter  windows  were  placed  in  the  combustor  so  that  the  combustion  process  could  be 
observed  and  the  steady  combustion  characteristics  established. 

The  torus  combustor  was  operated  over  the  following  range  of  conditions: 


Steady-state  chamber  pressure,  P  ,  Ib/sq  in.  abs . 40  to  120 

Oxidizer-flow  rate,  Ib/sec . 0.95  to  3.  5 

Fuel-flow  rate,  Ib/sec .  0.  15  to  0.  45 
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Total  flow  rate,  Ib/sec . 

Oxidizer -fuel  weight  ratio  . 

Oxidizer  injection  velocity,  ft/sec 

Hydrogen  injection  velocity,  v^  .,  ft/sec 

Oxidizer  pressure  drop,  AP’  . 

o 

Hydrogen  temperature,  F . 

Liquid-oxygen  temperature,  . 


1.  2  to  4.  0 

2.  0  to  8.  0 
15  to  55 

300  to  600 

0.  2  to  0.  8 
32  to  80 
.  .  -320 


To  excite  instability  a  tangential  blast  tube  10  feet  long  with  a  2 -inch  inside  diameter 
was  added  to  the  combustor.  A  burst  disk  that  ruptured  at  1000  pounds  per  square  inch 
was  placed  at  the  combustor  end  of  the  blast  tube.  The  tube  was  loaded  with  75  grams  of 
ammonium  perchlorate  containing  25  percent  by  weight  polybutadiene  acrylic  acid,  which 
was  ignited  during  the  run  and  which  provided  sufficient  gas  to  burst  the  diaphragm  and 
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Figure  10.  -  Experimental  and  theoretical  stability  boundaries  for  torus  combustor. 
Burning-rate  parameter,  0,015. 


produce  a  peak-to-peak  pressure  oscillation  of  20  percent  of  the  chamber  pressure  in  the 
torus  combustor. 

The  experimental  results  are  shown  in  figure  10,  where  the  disturbance  amplitude  is 

plotted  against  the  mass  accumulation  factor  27rM^/co.  ^t  .  For  the  unlike -triplet  and 

o  ^  w 

impinging-jet  injectors  the  20-percent  disturbance  always  damped.  The  instability  of  all 
runs  with  the  two  parallel-jet  injectors  indicated  that  the  combustion  noise  produced  dis¬ 
turbances  large  enough  to  excite  instability.  The  noise  level  was  assumed  to  be  the  same 
as  that  observed  for  the  triplet  and  impinging-jet  injectors  during  steady-state  combus¬ 
tion  (4  percent  of  chamber  pressure),  although  larger  disturbances  were  present  during 
ignition  and  while  flows  were  adjusting  to  their  average  rate.  The  instability  was  always 
initiated  sometime  after  ignition  but  before  average  flow  rates  were  achieved. 

The  curves  in  figure  10  are  cross  plots  of  data  presented  in  figures  4  and  6  for  a 
burning  rate  parameter  of  0.  015.  This  value  of  ^  is  for  the  torus  with  a  radius  PJi 
of  48  inches,  an  average  burning  rate  per  total  length  of  combustor  ni  of  0.  131,  and  a 
combustor  contraction  ratio  of  409  (obtained  by  dividing  the  volume  of  the  combustor 
by  the  product  of  the  nozzle  area  and  the  distance  from  the  injector  to  the  nozzle).  The 
theoretical  curves  show  that  the  vaporization  model  has  a  slightly  lower  instability  limit 
than  the  atomization  model  for  the  operating  conditions  of  the  torus  combustor.  The  ex- 
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perimental  value  of  the  mass  accumulation  factor  was  obtained  by  calculating  the  amount 
of  vaporized  propellant  as  a  function  of  length  (ref.  5).  These  calculations  agreed  with 
the  observation  at  the  window  that  all  of  the  propellant  was  vaporized  for  the  impinging- 
jet  and  triplet  injectors,  while  only  a  small  fraction  was  vaporized  with  the  parallel-jet 
injectors.  Thus  the  parallel-jet-type  injectors  provided  sufficient  liquid  in  the  combus¬ 
tor  to  significantly  reduce  the  calculated  disturbance  amplitude  required  to  excite  in¬ 
stability. 

The  experimental  and  calculated  resvilts  shown  in  figure  10  indicate  that  good  agree¬ 
ment  exists  between  theory  and  experiment.  Theory  indicates  that  a  disturbance  greater 
than  40  percent  is  required  to  excite  instability  with  the  impinging-jet  and  triplet  injec¬ 
tors,  and  experiments  showed  that  a  20-percent  disturbance  would  decay.  Similarly,  the 
parallel-jet-type  injectors,  which  were  always  unstable,  should  require  only  a  2-  to 
6-percent  disturbance  to  excite  instability.  Operating  conditions  of  the  torus  combustor 
did  not  completely  agree  with  the  conditions  used  in  the  calculations.  However,  these 
would  not  significantly  change  the  calculated  curves  because  the  calculations  are  quite 
insensitive  to  /,  v’ ,  AP',  and  y.  The  minimum  disturbance  is  very  sensitive  to  Av’; 
however,  the  value  of  0.  01  for  Av’  is  very  reasonable  in  that  it  represents  the  turbu¬ 
lence  level  available  as  a  driving  source  when  gas  and  liquid  velocities  are  very  low. 

High-speed  motion  pictures  taken  of  the  sprays  during  instability  indicated  that  the 
sprays  were  undergoing  breakup  due  to  crossflow  at  the  high-pressure  portion  of  the 
oscillations.  It  was  impossible  to  determine  the  rate  of  breakup  or  vaporization,  how¬ 
ever,  because  the  liquid-oxygen  sprays  are  not  clearly  outlined  in  the  pictures;  there¬ 
fore,  the  volume  of  liquid  could  not  be  measured  accurately  enough  to  obtain  slopes  and 
hence  atomization  or  vaporization  rates.  The  amount  of  liquid  in  the  combustor  varied 
during  the  cycle  as  predicted  by  the  theory. 


SUMMARY  OF  RESULTS 

Theoretical  calculations  have  shown  that  the  physical  processes  of  vaporization  and 
atomization  probably  are  the  controlling  parts  of  the  combustion  process  in  determining 
combustion  instability.  Chemical  reaction  is  not  likely  to  be  important  for  most  rocket 
combustors  because  there  is  insufficient  mass  or  energy  available  in  gaseous  form  to 
sustain  an  oscillation  even  though  the  reaction  rate  is. very  sensitive  to  pressure,  veloc¬ 
ity,  or  temperature  variations. 

Theory  also  indicates  that  the  stability  of  a  combustor  is  very  dependent  on  the  mass 
of  liquid  in  the  combustor  and  the  existence  of  any  vortex  velocity,  as  well  as  the  veloc¬ 
ity  difference  and  the  burning  rate  parameter  (ref.  2). 


17 


Limited  experimental  studies  with  a  combustor  similar  to  the  one  used  in  the  theory 
have  confirmed  the  stability  regions.  Experiments  also  confirm  qualitatively  the  vari¬ 
ations  that  occur  during  an  oscillation. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  February  9,  1965. 


18 


APPENDIX  A 


SYMBOLS 


cross-sectional  area  of  com¬ 
bustor,  sq  in, 

average  surface  area  of  drops, 
sq  in. 

injection  hole  area,  sq  in. 

nozzle -throat  area,  sq  in. 

combustor  contraction  ratio, 
Ac/At,  dimensionless 

average  speed  of  sound  in 
gases,  in.  /sec 

coefficient  of  discharge  of  in¬ 
jector,  dimensionless 

average  steady- state  mass 
fraction  of  unburned  gases, 
dimensionless 

characteristic  exhaust  veloc¬ 
ity,  ft/sec 

molecular  diffusion  coeffi¬ 
cient,  sq  in.  /sec 

activation  energy,  in.  Ib/lb 
mole 


function  of  gamma, 

9  \r+i/r-i 


r  +  1, 


acceleration  due  to  gravity, 
386.  09  in.  /sec^ 


viscous -dissipation  parameter, 
* 

/KF 

,0‘ 


7RP, 


k 

preexponential  constant  in 
Arrhenius  equation, 
cu  in.  /(lb  mole) (sec) 

fjP 

burning-rate  parameter, 

Rm/ji 

M 

mass  concentration,  Ib/cu  in. 

mass  concentration  of  unatom¬ 
ized  propellant,  Ib/cu  in. 

mass  concentration  of  unre¬ 
acted  propellant,  Ib/cu  in. 

mass  concentration  of  un vapor¬ 
ized  propellant,  Ib/cu  in. 

m 

average  burning  rate  of  propel¬ 
lant  per  length  of  combustor, 
fraction/in. 

molecular  weight  of  propellant, 
lb  mass/lb  mole 

Nd 

number  of  liquid  drops  per  unit 
volume  in  combustor, 
l/cu  in. 

chamber  pressure,  Ib/sq  in. 

O 

o 

steady-state  chamber  pressure 

p 

c,  max 

maximum  chamber  pressure 

P 

c,  mm 

minimum  chamber  pressure 

nondimensional  peak-to-peak 
pressure  disturbance, 

P  -  P  . 

c,  max  c,min 

^C,  0 

Pt 

tank  pressure,  Ib/sq  in. 

axial  gas  yelocity,  in.  /sec 

Py 

vapor  pressure  of  liquid, 

tangential  gas  yelocity,  in.  /sec 

Ib/sq  in. 

We 

Weber  number,  dimensionless 

AP’ 

nondimensional  steady-state 

z 

axial  coordinate 

pressure  drop  across  injector, 

Pp  ■  Pc,  o 

^C,  0 

z' 

a 

nondimensional  axial  coordinate, 
z/R 

correction  factor  for  mass  trans¬ 

R 

radius  of  combustor,  in. 

fer,  dimensionless 

^/i 

universal  gas  constant,  18  510 

7 

specific-heat  ratio,  dimensionless 

Re 

(in.  )(lb)/(^R)(lb  mole) 

Reynolds  number,  dimension¬ 
less 

e 

angular  coordinate 

gas  yiscosity,  lb/(in.  )(sec) 

average  drop  radius,  in. 

p 

gas  density,  Ib/cu  in. 

Sc 

average  jet  radius,  in. 

CO 

burning  rate,  lb/(sec)(cu  in. ) 

Schmidt  number,  dimensionless 

atomization  rate,  lb/(sec)(cu  in. ) 

T 

temperature,  ^R 

injection  rate,  lb/(sec)(cu  in. ) 

t 

time,  sec 

chemical  reaction  rate. 

t* 

breakup  time  of  liquid,  sec 

average  wave  time,  27rR/a,  sec 

nondimensional  time,  t/t 

lb/(sec)(cu  in. ) 

cOy  vaporization  rate,  lb/(sec)(cu  in. ) 

Subscripts: 

U 

w 

velocity  difference  between  liquid 
and  gas,  in.  /sec 

f 

£ 

fuel 

liquid 

V 

gas  velocity,  in.  /sec 

0 

steady-state  conditions 

Av 

steady-state  average  axial  ve¬ 

ox 

oxidizer 

locity  difference  between 
liquid  and  gas,  in.  /sec 

Superscript; 

v‘ 

nondimensional  gas  velocity,  v/a 

T 

nondimensionalized  parameter. 

Av‘ 

nondimensional  axial  velocity 

equal  to  parameter  diyided  by 

difference,  Av/a 

steady-state  average  vortex  ve¬ 

steady-state  parameter,  i.  e. , 

T'  =  T/Tq  (except  for  yelocities, 
which  are  diyided  by  speed  of 

v’ 

vr 

locity,  in.  /sec 

nondimensional  vortex  velocity. 

sound,  and  time,  which  is 
diyided  by  waye  time) 

Vyr/^ 
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APPENDIX  B 


EQUATIONS  FOR  RATES  AND  MASSES 
Injection  Rate 

The  local  instantaneous  injection  rate  is  given  by 

=  Cj^A.  |/(P^  -  P^)2p^g  (Bl) 

The  steady- state  injection  rate  is 


Atomization  Rate 


The  time  to  atomize  a  jet  or  drop  of  liquid  in  a  cross  current  for  values  of 
We/Re^'^^  much  larger  than  1  is  (ref.  6) 


t 


B 


".i 

^  (2We)^/^ 


(B4) 


where  Re  =  rj(Up/fJ.)  and  We  =  rj(pU  /a).  If  the  instantaneous  atomization  is  the  same 
as  the  average  for  complete  brea!^p,  the  atomization  rate  is 


a 


and 
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iL  l-25 


5/12 

=  M'(pT^^  — - 

Si  9 

(Av)"^ 


0.  625 


M'(pt)Vl2  j 

ci 


0.  625 


Vaporization  Rate 


The  rates  of  vaporization  for  large  Reynolds  numbers  as  given  in  reference  1  are  as 


follows: 


N  ,Dyy/,A jCkP 

a,  =_d - 1.  d_._.v  0  6 

/y/T2r^ 


1/2  roT- 


K<''  -  ''j.o)] 


N,D;;^,A.  ^Q!P„  1 /q 

=  d _ f  d,o _ Vr.  fiopl/J 


v,o  ;y/T2r 


0.6ScV’^[2r^^^(v  .-  v^^^) 


1/2/p 


-V  _  ^d  r.  o  -  /p 


^  ‘%,o  ^d,oVrd 


1/2  1/2 


'"o  ■  "^1,0 
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..211/4 


(BIO) 


Chemical  Reaction  Rate 

The  chemical  reaction  rate  is  given  by  a  second-order  rate  expression  in  which  an 
Arrhenius  rate  constant  is  assumed; 


u)  M  M  f 
r  _  r,  ox  r,  f 


‘^ox.f  “ox  “f 


k  exp 


(Bll) 


w  .  M„  „„  „  M, 


iVX  J-vx  r  ^ 

j_2.  -  r,ox,  o  r,f,  O]. 


m  X  iri 
ox, f  ox 


(B12) 


_ ^ _ =  M*  M'  „  exD  — ? —  /l  -  — 


=  M*  M'  f  exp  -- —  1  -  -j 
r,ox  r,i  y  -p 


(B13) 


Mass  of  Unatomized  Propellant 

The  mass  of  unatomized  propellant  is  determined  by  a  mass  balance; 


M  =  ^  +  /*  (o).  -  oj  )dt 

a  a,  o  1  a 


M’  = 


=  1  + 


^a,  o  X  \^i,  o  '^1,0/ 


t  dt' 
w 


(B14) 


(B15) 


Assuming  that  all  the  propellant  is  consumed  in  the  combustor  gives 


CO  =  CO.  =  OJ  ^ 

v,o  1,0  a,  o 
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a.  o 


(B16) 


Mass  of  Unvaporized  Propellant 


The  mass  of  unvaporized  propellant  is  also  determined  by  a  mass  balance; 


M  =  M 

V  V 


,0+/ 


(B17) 


M 

M’  = — ^ 
“v.O 


).  /*  /o)  W  \ 

/  L-J_  -^\ 

I  /  Vw.  O).  / 

V,  o  ^0  \  1, 0  1, 0/ 


t  dt' 
w 


(B18) 


CO  ~  U).  =  CO 

V,  0  1,  o  a,  o 


(B19) 


m;  =  1  + 


rV 

t  /  (co'  -  coMdt' 

w  .A  '  a 


(B20) 


Mass  of  Unreacted  Propellant 


The  amount  of  unreacted  propellant  is  also  determined  by  a  mass  balance; 


M  =  M  ^  +  A  (co  -  CO  )dt 
r  r,o  L  ^  V  r^ 


{B21) 


where 


M  -  P 
r, 0  ’  0^0 


(B22) 
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(B23) 


since 


O).  t  CO.  ^271’R 

=  -iiO - =  277,'?f(y) 


(See  eqs.  (A29)  and  (A31)  of  ref.  2. ) 


m;  =  1  + 


2-!i^i{y) 


/' 

•'A 


(co'  - 

'  V 


copdt' 


{B24) 
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